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SUMMARY 

I  An  analytical  study  of  a  flat  condenser -radiator  with  a  double  fin-tuhe 
geometry  (closed  sandwich)  with  variable  tube  side-wall  thickness  was  performed 
for  a  Rankine  space-power  electric-generating  system.  The  analysis  of  the 
double  fin  radiator  included  consideration  of  tube  and  header  pressure  drops, 
meteoroid  protection  for  the  tubes  and  headers,  along  with  a  detailed  presenta¬ 
tion  of  the  heat  rejection  analysis  and  total  weight  characteristics.  The 
double  fin-tube  radiator  is  compared  to  a  conventional  central  fin-tube  con¬ 
figuration  on  a  heat  rejection  to  weight  basis  for  a  four-panel  radiator  con¬ 
figuration. 


Both  fin  and  tube  geometries  are  compared  on  the  basis  of  the  same  power 
level,  working  fluid  temperature,  tube  and  header  pressure  drop,  radiator  mate¬ 
rial,  and  meteoroid  protection  criteria.  A  beryj^/um  radiator  for  a  1-megawatt 
"^system  and  a  coluijijpium  alloy  radiator  for  a  500-'fcflowatt  system,  both  at  a  ra^;; — 
a  rating  temuerat^  of  1700°  R,  were  chosen  for  the  weight  and  geometry  com¬ 
parisons  .\  _  ^ 


The  conclusion  reached  indicates  a  substantial  weight  savings  can  be  real' 


ized  with  the  double  fin-tube  arrangement  if  the  tube  side-wall  thickness  can 
be  reduced  as  a  result  of  a  possible  meteoroid  bumper  effect  of  the  enclosing 
fins.  Weight  reductions  compared  to  the  central  fin-tube  geometry  of  up  to 
32  to  39  percent  were  shown  to  be  possible  for  the  maximum  reduction  in  side- 
wall  thickness  in  the  two  examples  considered.  This  result  fixrther  substanti¬ 
ated  the  preliminary  conclusions  given  in  an  earlier  reference  that  compared 
the  double  and  central  fin-tube  configurations  neglecting  the  effects  of 
headers,  pressure  drops,  tube  wall  temperature  drop,  and  powerplant  thermo¬ 
dynamic  cycle  considerations.  Thus  an  incentive  is  offered  to  further  investi¬ 
gate  the  meteoroid  b\uiiper  screen  concept  and  its  application  in  the  double  fin- 


tube  geometry  radiator. 


INTRODUCTION 

The  generation  of  large  amounts  of  electric  power  in  space  using  the 
Rankine  cycle  concept  requires  that  a  large  amount  of  waste  heat  be  rejected 


from  the  vorking  fluid.  This  vaste  heat  is  the  amount  of  energy  that  must  he 
rejected  from  the  working  fluid  vapor  leaving  the  turbine  in  order  to  complete¬ 
ly  condense  it.  Since  radiation  is  the  only  mode  of  heat  transfer  for  reject¬ 
ing  this  energy  and  since  maximum  temperatures  are  limited^  the  resulting  ra¬ 
diator  surface  areas  and  weights  are  generally  large.  Previous  studies  on 
space  radiators  (refs.  1  to  5)  have  indicated  that  a  central  fin-tube  arrange¬ 
ment  is  feasible  since  it  reduces  radiator  weight  by  reducing  the  prime  surface 
area  vulnerable  to  critical  damage  by  impacting  meteoroids.  Reference  6  pre¬ 
sents  the  analysis  and  results  of  a  central  fin-tube  radiator  using  flat  plate 
fins  of  constant  cross  section. 

A  preliminary  comparison  of  several  fin-tube  configurations  was  carried 
out  in  reference  7  that  did  not  include  cycle  considerations^  vapor  and  liquid 
headers^  and  tube  and  header  pressure  drops.  The  results  of  reference  7  indi¬ 
cated  that  a  substantial  weight  savings  could  be  obtained  by  using  a  double  fin 
tube  with  reduced  rectangular  tube  side-wall  thickness  instead  of  the  central 
or  open  sandwich  fin-tube  arrangements.  A  reduced  tube  side-wall  thickness  can 
result  ^ from  the  assumption  that  the  two  fins  of  the  double  fin-tube  configura¬ 
tion  will  act  as  a  bumper  screen  to  incident  meteoroids^  thus  allowing  some 
reduction  in  the  thickness  of  the  tube  side  wall.  The  double  fin-tube  configu¬ 
ration  is  attractive  from  a  structural  viewpoint  since  it  provides  a  rigid 
structure  that  has  a  continuous  smooth  surface  that  could  be  the  vehicle  skin. 
These  attractive  features  motivated  the  need  for  a  more  sophisticated  analysis. 

The  purpose  of  this  study  is  to  analyze  in  a  more  comprehensive  manner  the 
heat  rejection  and  weight  characteristics  of  a  double  fin-tube  configuration 
with  variable  rectangular  tube  side-wall  thickness^  and  to  more  accurately 
identify  the  potential  weight  savings  over  radiators  with  a  central  fin-tube 
geometry.  The  vapor  and  liquid  headers  are  also  included  in  the  geometry  ana¬ 
lyzed  since  they  directly  affect  the  meteoroid  protection  req.uirements  of  the 
radiator  tubes  and  can  be  a  significant  portion  of  the  total  radiator  weight. 
There  is  also  a  sizable  amount  of  heat  that  can  be  rejected  from  the  vapor 
header  for  large  power  systems  (ref.  l).  Working  fluid  pressure  drops  were 
also  considered  for  the  tubes  and  headers. 

A  rectangular  cross-section  fin  was  chosen  for  the  analysis  and  compari¬ 
son.  A  one -dimensional  approach  was  taken  in  the  development  of  the  fin  energy 
balance  equation  with  the  assimiption  that  the  base  temperature  of  the  fin  is 
equal  to  the  outer  surface  temperature  of  the  tube  armor.  The  radiator  heat 
rejection  and  weight  analysis  was  carried  out  for  two  typical  Rankine  power- 
planets;  a  1-megawatt  potassium  cycle  with  a  1700*^  R  beryllium  aimor  and  fin 
radiator^  and  a  500-kilowatt  potassium  cycle  with  a  1700^  R  columbium  alloy 
radiator.  This  report  presents  the  results  of  the  heat  rejection  and  weight 
analysis  for  the  double  fin-tube  radiator  along  with  a  comparison  of  these  re¬ 
sults  with  those  given  in  reference  6  for  the  central  fin-tube  configuration. 
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radiator  planform  area^  sq  ft 
vulnerable  area,  sq  ft 
penetration  correction  factor 
constant 

specific  heat,  Btu/(lh)(^F) 

sonic  velocity  in  armor  material,  V^aS/Pa^  sec 
tube  diameter,  ft 

young  ^s  modulus  of  armor  material,  Ib/sq  ft 

angle  factor,  fraction  of  the  energy  leaving  a  surface  that  is  incident 
upon  another  surface 

vapor  header  occlusion  factor 

units  conversion  factor,  32.17  ft/sec^ 

incident  energy,  Btu/(hr)(sq  ft) 

heat  of  condensation,  Btu/lb 

mechanical  equivalent  of  heat,  778  (ft)(lb)/Btu 
factor,  (q-jj  +  q^) 

fluid  turning  loss  factor  from  header  to  tubes 

thermal  conductivity,  Btu/(hr )  (ft )  (^R) 

minim-um  half  length  of  fin,  equal  to  L*  -  Rq,  ft 

one-half  the  tube  center  to  center  distance,  ft 

actual  half  length  of  fin,  ft 

number  of  radiator  tubes 

conductance  parameter,  aZ^T-^/kt 

cycle  fluid  pressure,  Ib/sq  ft 

poverplant  output,  kw 

probability  of  zero  punctures 

heat  rejection  rate,  Btu/hr 


3 


% 

% 

*^ej 

Qualj 

q. 

R 

% 

£ 

Re 

T 

IP* 

t 


tube  radiant  heat  rejection  rate,  Btu/hr 

fin  radiant  heat  rejection  rate  for  a  fin  length  2Z 
both  sides,  Btu/hr 

total  radiator  heat  rejection  rate,  Btu/hr 

vapor  header  entrance  quality 

vapor  header  heat  rejection  rate,  Btu/hr 

heat  rej.ection  per  unit  length  of  tube,  Btu/hr  ft 

radius ,  ft 

tube  side  wall  to  tube  centerline  dimension,  R-j^  =  R^ 

fraction  of  flow  area  occupied  by  one  phase 
Reynolds  number 
temperature ,  °R 

vapor  saturation  temperature  at  tube  inlet,  °R 
thickness  of  fin,  ft 


radiating  from 


u  velocity  of  vapor,  ft /sec 

V  velocity  of  liquid,  ft/sec 

V  average  meteoroid  velocity,  ft/sec 

W  weight,  lb 

W  weight  flow  per  tube,  lb /sec 

W  panel  width,  ft 

X  normalized  distance  coordinate,  x/z 

fraction  of  total  heat  rejected  by  the  vapor  header 
Xj.^  fraction  of  total  heat  rejected  by  the  tubes  and  fins 
X  coordinate  measuring  distance  along  lower  fin,  ft 

V  normalized  distance  coordinate,  y/z 

y  coordinate  measuring  distance  along  upper  fin,  ft 
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z 


radiator  tvibe  lengthy  ft 


constants  in  penetration  formula 
&  tute  -wall  thickness^  ft 

6  surface  hemispherical  emissivity 

Tj*  thermal  effectiveness 

e  normalized  temperature^  T/T-^ 

|j.  viscosity^  Ih/ft-sec 

p  density;,  Ih/cu  ft  (unless  otherwise  specified) 

c  Stefan-Boltzmann  constant^  0.173x10“®  Btu/(hr)(sq 

T  mission  exposure  time^  days 

two-phase-flow  parameters 
Subscripts : 
a  armor 

h  tube  base  sixrface 

c  tube  liner 

cond  conduction 
F  friction 

f  fin 

g  vapor  phase 

i  inside 

LH  liquid  header 

£  liquid 

m  momentum 

o  outside 

P  particle 


R 


total 


rad  radiation 

s  side  -wall 

t  tulDe 

tot  total  flow,  liquid  and  vapor 
VH  vapor  header 

X  normalized  distance  coordinate,  x/z 

X  coordinate  measuring  distance  along  lower  fin,  ft 
Y  normalized  distance  coordinate,  y/z 

y  coordinate  measuring  distance  along  upper  fin,  ft 
0  conditions  at  tube  inlet 

1  hase  surface  1 

2  hase  surface  2 

3  radiator  inlet  conditions 


RADIATOR  COITFIGURATION  AM)  THERIVIODYNAraC  CYCLE 

The  general  radiator  panel  configuration  considered  for  the  analysis  is 
shown  in  figure  1  where  the  vapor  from  the  turbine  exhaust  is  distributed  to 
the  finned  tubes  by  a  vapor  header.  The  heat  radiated  from  the  vapor  header 
and  finned  tubes  causes  the  vapor  to  condense.  The  condensate  is  then  sub¬ 
cooled  and  collected  in  the  liquid  header  before  being  sent  to  the  condensate 
pump.  This  scheme  can  be  modified  by  dividing  the  radiator  into  a  number  of 
nonredundant  segments,  each  of  which  could  be  treated  as  a  separate  entity. 

The  detailed  cross-section  drawing  of  the  double  fin-tube  with  variable 
tube  side  wall  composed  of  tube  liner  inserted  in  an  armor  block,  which  pro¬ 
vides  meteoroid  protection,  and  two  rectangular  fins  is  shown  in  figure  2.  The 
tube  liner,  which  is  exposed  to  the  working  fluid,  must  be  capable  of  with¬ 
standing  possible  corrosion.  The  liner  thickness  must  also  be  compatible  with 
current  fabrication  capabilities  and  structural  requirements,  and  the  presence 
of  the  liner  may  also  substantially  reduce  the  required  armor  thickness 
(ref.  8).  The  liner  thickness  is  arbitrarily  taken  as  8^,  =  0.040  D.  with 
a  minimum  wall  thickness  set  at  0.020  inch.  The  liner  thickness  was  increased 
as  the  inside  diameter  was  increased,  so  as  to  provide  necessary  stiffening  and 
strengthening  of  the  radiator  tubes. 

For  the  double  fin  geometry  of  figure  2,  the  liner  can  be  damaged  by  im- 
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pacting  meteoroids  in  two  general  ways.  The  first  is  hy  any  primary  impacts 
occurring  on  the  outer  exposed  surfaces  of  the  tube  block.  These  impacts  are 
assumed  to  obey  the  conventional  armor  penetration  and  damage  relations  devel¬ 
oped  for  tubes  (refs.  8  to  10)^  with  vulnerable  area  given  by  4R^ZN.  Accord¬ 
ingly,  the  armor  thickness  5a  vas  determined  using  the  criterion  of  refer¬ 
ence  9.  The  armor  thickness,  which  is  a  result  of  the  optimization  program 
(ref.  6),  is  applied  in  full  on  the  upper  and  lower  surface  of  the  tube.  A 
second  damage  source  can  arise  from  a  spray  of  particles  on  the  armor  block 
side  surface  ((Dq  -  2t)ZN)  resulting  from  impacts  on  the  fin  surfaces.  In  view 
of  the  bumper  action  involved  and  the  obliquity  of  the  secondary  impacts,  how¬ 
ever,  a  reduction  will  undoubtedly  be  allowed  in  the  armor  thickness  required 
by  the  tube  block  side  wall  to  resist  the  effects  of  these  secondary  impacts. 
Since  no  specific  relations  are  at  present  available  for  the  determination  of 
this  side-wall  thickness,  a  parametric  variation  of  the  6g/5^  is  used  to  ex¬ 
amine  the  effects  of  reduced  side-wall  thickness  on  radiator  weight  and  geom¬ 
etry.  The  parameter  defined  as  the  fraction  of  the  armor  thickness 

retained  on  the  enclosed  side  of  the  tube  block. 

The  vapor  header  takes  the  form  of  a  hollow  paraboloid  whose  wall  consists 
of  0.12-inch-thick  liner,  which  was  arbitrarily  chosen,  and  whose  meteoroid 
armor  protection  thickness  is  the  same  as  that  required  by  the  tubes.  The  par¬ 
abolic  shape  insures  constant  velocity  in  the  header.  For  simplicity,  the 
liquid  header  was  designed  with  a  constant  diameter  and  a  fluid  velocity  of 
4  feet  per  second,  so  that  a  very  low  press-ure  drop  would  result.  The  liner 
for  the  liquid  header  follows  the  same  schedule  with  the  inside  diameter  as  do 
the  tube  liners.  However,  a  maximum  liner  thickness  is  set  at  0.12  inch.  The 
liquid  header  also  has  meteoroid  armor. 

The  assumptions  given  herein  for  the  double  fin-tube  radiator  geometry, 
with  the  exception  of  the  fin-tube  configuration,  also  hold  for  the  central  fin 
tube  that  was  analyzed  in  reference  7  and  shown  in  figure  3. 

The  thermodynamic  cycle  used  in  this  analysis  is  the  Rankine  cycle,  which 
uses  a  working  fluid  that  undergoes  a  change  of  phase.  The  working  fluid  is 
condensed  in  the  radiator,  which  results  in  a  near  isothermal  condition  pre¬ 
vailing  in  the  tubes  and  vapor  header.  In  order  to  show  sample  results  and 
compare  the  two  fin-tube  geometries  on  a  heat  rejection  per  unit  weight  and  on 
a  radiator  geometry  basis,  two  power  levels  were  chosen.  Potassium  was  chosen 
as  the  working  fluid  in  the  cycle  for  both  power  levels  with  a  peak  turbine 
inlet  temperature  of  2460^  R  and  a  radiator  temperature  of  1700^  R.  It  was 
also  specified  that  the  radiator  tubes  would  subcool  the  working  fluid  100*^  R. 
Additional  cycle  requirements  such  as  turbine  and  generator  efficiencies  were 
set  at  0.75  and  0.90,  respectively,  with  10  percent  of  the  generator  output  re¬ 
quired  for  accessories  and  controls.  The  emittance  of  the  radiator  was  taken 
to  be  0.90,  and  the  effective  sink  temperature  for  the  radiator  was  assimied  to 
be  0^  R.  The  foregoing  values  plus  the  cycle  temperatures,  working  fluid,  and 
the  power  level  chosen  for  the  analysis  and  comparison  enabled  the  determina¬ 
tion  of  the  total  heat  rejection  rate  and  the  mass  flow  rate  of  the  working 
fluid  for  the  radiator  design  inputs.  Additional  information  required  from  the 
cycle  analysis  is  the  quality  of  the  working  fluid  entering  the  vapor  header. 

The  analysis  of  the  thermodynamic  cycle  used  is  given  in  detail  in  reference  6. 
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HEAT-TRMSZER  ANALYSIS 


Approach  and  Assixaiptions 

The  analysis  considers  the  general  case  of  two  rectangular  profile  fins  of 
length  I  attached  to  a  tube  enclosed  in  an  armor  "block  forming  a  double  fin 
configuration  as  shown  in  figure  2.  Energy  input  to  the  fin  is  comprised  of 
heat  conduction  along  the  fin  from  the  two  tube  side-wall  surfaces.  For  any 
specific  choice  of  L/R^  the  fin  length  I  as  shown  in  figure  2  will  depend 
on  the  value  of  tube  side-wall  thickness. 

The  specific  assumptions  used  in  the  development  of  the  heat-transfer  re¬ 
lations  for  the  double  fin-tube  geometry  as  well  as  the  central  fin-tube  geome¬ 
try  are 

(1)  Radiator  emittance  is  constant  with  temperature. 

(2)  For  the  determination  of  radiator  temperature  variation  the  radiator 
surfaces  act  as  blackbodies  with  incident  and  emitted  radiation  governed  by 
Lambert ^s  cosine  law. 

(3)  Hemispherical  radiation  to  space  from  both  outer  surfaces  of  the  radi¬ 
ator  to  a  0^  R  space  sink  temperature. 

(4)  The  base  surface  temperature  of  the  fin  is  assumed  constant  along  the 
length  of  the  tube  and  eq_ual  to  the  tube  block  outer  surface  temperature. 

(5)  Steady-state  one -dimensional  heat  flow  exists  in  the  fins  and  tube 
block. 

(6)  Material  properties  are  constant  along  the  length  of  fin  and  tube 
block  and  are  evaluated  at  the  fin  base  temperature. 

(?)  The  development  of  the  fin  and  tube  angle  factors  for  the  radiant  in¬ 
terchange  between  fin  and  tube  side  walls  is  based  on  an  infinite  length  of 
tube  and  fin^,  with  fin  thickness  assumed  negligible  in  the  development. 

(8)  The  inside  tube  wall  temperature  is  circumferentially  uniform  and 
eq_ual  to  the  stagnation  temperature  of  the  fluid  at  the  inlet  of  the  header. 


Tube  Wall  Temperature  Drop 


Since  the  tube  block  will  have  thick  walls  due  to  the  armor  required  for 
meteoroid  protection^  a  significant  temperature  drop  will  occur  across  the 
block.  The  assumed  heat -transfer  paths  and  the  various  fin  and  tube  block  tem¬ 
peratures  are  shown  in  figure  4.  The  inside  tube  wall  temperature  T*  is  de¬ 
termined  from  reference  6  as 


T*  = 


T.  1  - 


Jgh 


(1) 
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vhere  is  the  radiator  fluid  stagnation  temperature  at  the  header  inlet  in 

degrees  Rankine  and  Uq  is  the  tube  inlet  vapor  velocity.  The  turning  loss 
factor  in  equation  (l)  vas  set  at  1.15  (ref.  6). 

The  relation  hetveen  the  tube  inside  surface  temperature  and  the  tube 
block  outer  surface  temperature  is  based  on  a  simplified  approach  that  assumed 
heat  is  transferred  from  the  tube  inner  surface  to  the  exposed  surface  of  the 
tube  block  by  one -dimensional  conduction.  This  heat  path  vas  chosen  since  it 
represents  the  greatest  flow  of  energy.  The  heat  transmission  by  conduction  is 
assumed  to  travel  a  distance  6^^  through  a  cross-sectional  area  2E-^  dZ.  The 
expression  for  the  heat  conduction  is 

2kR. 

I-Scona  =—<“*-  (2) 

^a 

For  simplicity  It  is  assumed  that  the  temperature  drop  T*  -  is  gov¬ 
erned  primarily  by  the  radiant  heat  transfer  from  the  exposed  sOTface  of  the 
block  (neglecting  the  conduction  to  the  fins).  The  expression  for  radiation 
from  one  surface  element  may  be  vrltten  as 

■iW  =  (=> 

When  equations  (2)  and  (3)  are  combined  for  dq^^^^^  =  dq^^d^  resultant 
equation  for  the  approximate  temperature  drop  through  the  tube  block  wall  is 
given  as 

ae&a^b  +  k(Tb  -  T*)  =0  (4) 

The  temperature  T-j^  obtained  from  equation  (4)  is  then  used  in  the  de¬ 
velopment  of  the  fin  heat-transfer  relations.  The  temperature  of  the  tube  side 
vail  is  also  assumed  to  be  equal  to  the  temperature  of  the  tube  block  exposed 
surface  Ty^  regardless  of  the  thickness  of  the  tube  side  wall  (all  values  of 

^s/^a)- 


Fin  Temperature  Profile 

Formulation  of  eq-uations.  -  Considering  an  element  of  the  fin  surface  in 
figure  2  and  employing  the  previous  assumptions^  the  law  of  energy  conservation 
for  an  element  can  be  expressed  as  the  energy  balance  between  the  net  heat 
transfer  due  to  conduction  and  radiation  (ref.  ll) : 

^*^cond  ^^ad 

The  net  internal  heat  conduction  through  the 
length  Z  is  expressed  as 

^^cond  “  ^  ^-ktZ 


=  0 


(5) 


element  of  thickness  t  and 


dx> 


dx 


(6) 
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The  net  radiant  heat  rejection  from  both  sides  of  the  element  of  the  fin  is 
composed  of  its  emission  minus  the  incident  energy.  This  exi)ression  is  given 


'^^rad  ~  ( 


2oT^  -  dA^ 


The  incident  energy  term  in  equation  (?)  is  composed  of  the  incident  en- 

ergy  from  the  two  side  walls  to  an  element  of  area  on  the  fin 

oT,'^^,.  .  +P..  1  fo, 


hfdA^-1  +  ^dA^-2 


and  the  energy  leaving  an  opposing  fin  surface  incident  upon  the  other 


oT"^  dP , .  ^  . 

y 


Substitution  of  expressions  (8a)  and  (8b)  into  equation  (7)  yields 


2aT^  -  oT^IF 


X  “  ‘^b\^dA^-l  +  ^dA^-2)  + 


^y  ^dAy-dAy 


Introducing  equations  (6)  and  (9)  into  equation  (5)  and  the  dimensionless  vari- 
ahles 

’=-1 

y.2 


c  kt 


yield  the  expression 


The  angle  factors  in  equation  (lO)  are  evaluated  using  a  relation  (ref. 

12,  eq.  (31-58))  that  applies  to  parallel  surfaces  of  infinite  length.  For  the 
configuration  of  figure  2  the  angle  factors  in  question  are  given  as 
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Introduction  of  equations  (ll)^  (12)^  and  (13)  and  the  symmetry  of  the 
temperature  distribution  about  x  =  ly  that  is^  T(x)  =  T(2Z  -  x)  ^  into  equa¬ 
tion  (10)  yield  the  expression 
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Computational  procedure .  -  It  was  necessary  to  use  numerical  techniques  to 
solve  equation  (14)  for  the  temperature  profile  of  the  fin  of  length  Z.  The 
numerical  solutions  were  carried  out  on  an  lEM  7094  electronic  digital  com¬ 
puter  . 

The  second  order  differential  equation  d^e/dX^  =  +  f(x)  was  solved 

for  the  6  profile  in  the  interval  X  =  0  to  X  =  1.  The  boundary  conditions 
are  9=1  at  X  =  0  and  d0/dX  =0  at  X  =  1^  where  B  is  a  constant  and 
f(X)  is  a  function  of  X.  Kalaha's  method  was  used  and  is  now  described. 

The  term  6^  was  approximated  by  a  linear  function  of  9,  namely  by  the 
first  two  terms  of  its  Taylor  series  expansion.  Central  differences  were  used 
to  obtain  an  e^^ression  for  d^e/dX^.  At  each  point  of  a  mesh  on  the  interval 
X  =  0  to  X  =  1^  the  equation  d^s/dX^  =  B0^  +  f(X)  was  thus  approximated  as  a 
linear  function  of  0.  The  set  of  resulting  tridiagonal  linear  equations  was 
reduced  to  two-diagonal  form,  and  these  were  then  solved  for  the  9  profile  by 
backward  substitution.  An  initial  0  profile  guess  was  used,  and  each  succes¬ 
sive  iteration  yielded  approximately  one  decimal  place  of  accuracy.  Solutions 
were  obtained  as  a  function  of  the  input  parameters  and  Z/Rq. 

Temperature  profile  results.  -  Each  solution  of  the  fin  energy  equation, 
which  is  Independent  of  power  and  temperature  level  or  tube  inside  diameter, 
provided  a  temperatxire  distribution  along  the  fin.  Results  are  plotted  in  fig¬ 
ure  5  as  a  function  of  position  X  on  the  fin  for  several  parametric  values  of 
Ng  and  Z/Rq.  It  is  seen  that  the  temperature  drop  along  the  fin  is  very  small 
when  the  conductance  parameter  is  small,  which  Indicates  a  low  thermal 

resistance  of  the  fin.  Additionally,  there  is  little  effect  of  Z/Rq  when 
is  small,  although  a  greater  temperature  drop  occurs  with  increasing  Z/Rq 
ratio.  This  holds  regardless  of  the  choice  of  the  tube  wall  thickness  ratio 
since  only  the  1/Rq  ratio  is  considered  in  the  development  of  equa¬ 
tion  (14). 


Radiator  Effectiveness 

Fin  heat  rejection.  -  After  the  temperature  distribution  and  the  slope  of 
the  temperature  distribution  curve  at  X  =  0  have  been  determined,  the  net 
heat  transferred  by  the  fin  can  be  calculated.  The  rate  of  heat  loss  from  the 
fin  outer  surface  at  any  point  x  on  the  fin  from  a  differential  area  NZ  dx 
is  aeZNTi  dx.  The  overall  rate  of  heat  loss  Q^.  from  the  pair  of  fins  that 
forms  the  double  fin  tube  is 

/2Z  /•! 

eaZNTl  dx  =  4CTZNZT-^e  /  e^(X)dX  (16) 

•'Q 

This  equation  is  evaluated  using  the  solutions  of  the  fin  energy  equa¬ 
tion  (14).  Comparison  of  this  fin  energy  rejection  to  the  total  heat  loss  from 
both  sides  of  an  isothermal  fin-tube  section  of  length  2(Z  +  Rq,)  can  be  ex¬ 
pressed  in  dimensionless  terms  as 
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(17) 


n*  - 
rif  - 


% 


r 


e^(x)dx 


4€aZNz(l  +  ^j'^b 


where 


and  6(X)  is  a  function  of  X  for  specific  values  of  L/Rq,  and  N^. 

Ec^uation  (1?)  can  be  defined  as  the  thermal  effectiveness  of  the  double  fin. 

The  results  obtained  from  the  solutions  of  eq^uation  (17)  for  fin  effec¬ 
tiveness  are  presented  for  an  example  case  in  figure  6  as  a  function  of  l/R^ 
for  parametric  values  of  and  In  the  aforementioned  equations  de¬ 

scribing  fin  effectiveness  (eqs.  (17)  and  (18))^  the  tube  vail  thickness 
must  be  known  in  order  to  obtain  solutions.  In  order  to  obtain  the  actual 
value  of  the  tube  inside  diameter  tube  liner  thickness  b^^  power  and 

temperature  level,  materials,  meteoroid  protection  criteria,  tube  and  header 
pressure  drop,  and  definition  of  radiator  tube  vulnerable  area  must  be  speci¬ 
fied.  Thus,  solutions  for  equation  (17)  require  complete  solutions  for  the 
entire  radiator.  Such  radiator  solutions  will  be  described  in  detail  in  the 
section  RADIATOR  WEIGHT  MD  GEOMETRY. 

Inspection  of  the  curves  in  figure  6  reveals  that  for  any  fixed  value  of 
L/Rq  and  effectivenss  decreases  with  increasing  values  of 

The  fin  effectiveness  is  also  seen  to  decrease  with  increasing  bg/ba  con¬ 
stant  Rq-  When  l/Rq  equals  zero  the  value  of  fin  effectiveness  does  not 

equal  zero  when  the  value  of  bg/ba  is  other  than  1.  This  is  due  to  an  in¬ 

crease  in  the  fin  length  caused  by  the  reduction  of  tube  side-wall  thickness 
for  any  choice  of  L/Rq  ratio.  Thus,  at  l/Rq  equal  to  zero,  a  fin  remains 
of  length  Z  =  ba  -  bg. 

Tube  heat  rejection.  -  The  net  heat  loss  from  the  external  surface  of  the 
tube  block  is  just  its  radiant  emission  since  there  is  no  incident  energy  from 

other  parts  of  the  system.  This  energy  rejection  when  compared  to  the  total 

heat  loss  from  both  surfaces  of  an  isothermal  fin-tube  section  of  length 
2(Z  +  R-j^)  can  be  defined  as  the  thermal  effectiveness  of  the  exposed  surface 
of  the  tube.  This  expression  is 


^b 


1  + 


Rn 


(19) 
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where  Z/r^  is  defined  in  terms  of  Rq^  L,  63,^  ^s/^a  eq.uation  (18), 

Th©  1163,13  PS/diaiiion  fpoiu  iiliB  IiuTdc  sido-'WB,!!  surfaces  is  already  included  as  a 
contribution  to  the  heat  radiated  from  the  fin^  since  in  effect  the  radiant 
interchange  acts  in  a  manner  shnilar  to  the  heat  conduction  dovn  the  fin.  This 
aspect  is  fully  explained  in  reference  11.  Solutions  for  equation  (19)  also 
require  values  of  armor  thickness  6g^  and  thus  are  a  result  of  the  radiator 
calculation  procedure  for  the  double  fin-tube  block  configuration. 

Results  of  equation  (l9)  for  tube  effectiveness  are  shown  plotted  for  an 
example  case  in  figure  7.  Tube  thermal  effectiveness  is  seen  to  decrease  as 
L/Rq  is  increased^  which  indicates  a  decrease  in  the  relative  importance  of 
the  tube  portion  as  the  fin  length  increases.  Decreasing  the  tube  wall  ratio 
^s/^a  ^Iso  results  in  a  reduced  tube  effectiveness^  which  is  caused  by  a  re¬ 
duction  in  the  tube  block  outer  surface  area  and  thus  the  energy  it  can  reject 
to  space.  It  is  also  seen  from  the  figure  that  for  an  L/Rq  ratio  equal  to 
zero  the  tube  effectiveness  does  not  equal  1  if  the  ratio  is  less  than  1. 

This  is  brought  about  because  a  fin  of  length  6^  -  5^  remains  (see  eq.  (15)). 

Total  fin-tube  heat  rejection.  -  For  calculation  purposes^  it  is  desirable 
to  formulate  the  total  fin-tube  thermal  effectiveness  that  can  be  determined  by 
summing  the  results  of  equations  (17)  and  (19)  for  fin  effectiveness  and  tube 
effectiveness^  respectively.  This  expression  is 


* 

'Qr  - 


Of  +  Qb 


(20) 


Figure  8  shows  a  plot  of  total  fin-tuhe  effectiveness  against  L/Rq  ratio 

for  several  choices  of  conductance  parameter  and  tube  block  side-wall 

ratio  Sg/Sa  example  case.  Inspection  of  the  curves  shown  in  figure  8 

reveals  that  for  any  fixed  value  of  L/Rq  and  the  total  fin-tuhe  thermal 

effectiveness  decreases  as  the  Sg/h^  ratio  decreases.  This  is  reasonable 
since  as  the  Sg/Sa  ^atio  decreases  the  amount  of  isothermal  base  surface  de¬ 
creases.  The  fin-tube  effectiveness  is  also  reduced  by  increasing  the  conduc¬ 
tance  parameter  or  by  increasing  the  L/Rq  ratio. 

Negligible  variation  in  the  magnitude  of  r^*  was  observed  for  the  two 
power  levels  chosen  for  the  comparison.  Thus  the  curves  given  in  figure  8  for 
the  500-kilowatt  case  also  apply  for  the  1-megawatt  power  level  case.  It  was 
also  noticed  that  variations  of  the  tube  inside  diameter  had  no  appreciable  ef¬ 
fect  on  the  value  of  T)g  at  a  specific  value  of  the  L/Rq  ratio,  the  Bg/ba 
ratio,  and  the  conductance  parameter 


PRESSURE  DROP  CONSIDERATIONS 

Another  factor  that  is  required  to  determine  the  geometry  and  weight  of  a 
radiator  is  the  pressure  drop  determination  in  the  radiator  tubes  and  headers. 
This  aspect  of  radiator  design  helps  to  determine  the  vapor  header  geometry  and 
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the  required  tube  diameter^  the  tube  lengthy  and  the  number  of  tubes.  The 
equations  presented  are  for  a  Rankine  cycle  condenser-radiator  vhich  assumed 
that  vapor  flov  in  the  vapor  header  to  be  of  the  same  quality  as  the  turbine 
exhaust^  tvo-phase  flov  in  the  radiator  tubes  and  all-liquid  flov  in  the  liquid 
header.  The  development  of  the  equations  given  in  this  section  are  given  in 
reference  6 . 


Vapor  header.  -  The  determination  of  the  pressure  drop  in  the  vapor  header 
is  simplified  by  assuming  that  only  the  gas  phase  affects  the  pressure  drop. 
This  pressure  drop  is  expressed  as  a  ratio  of  to  the  header  inlet  pres¬ 

sure  P3  "With  the  resulting  ratio  kept  constant  for  comparative  purposes. 

This  equation^  which  is  for  turbulent  flow^  is 


0.00357  PgU^gW 


(21) 


where  Uy;^  is  the  uniform  vapor  velocity  in  the  parabolic  header  (based  on  the 
turbine  exhaust  quality  and  neglecting  the  flow  area  occupied  by  the  liquid) 
and  Re  is  the  vapor  Reynolds  number  based  on  the  vapor  header  maximum  diam¬ 
eter  Pyij-  The  value  of  is  obtained  from  the  expression 


2W  Qual^ 


1/2 


(22) 


and  the  total  panel  width  W  in  equation  (21)  Is 


(23) 


Odie  term  Qual3  In  equation  (22)  Is  the  vapor  quality  at  the  entrance  to  the 
vapor  header.  The  number  of  condensing  tubes  R  In  the  previous  equation  is 
determined  from  the  tube  pressure  drop  analysis  in  conjunction  with  the  opti¬ 
mization  procedure  used  for  this  fin-tube  geometry. 

The  amo\mt  of  heat  rejected  to  space  from  the  parabolic  vapor  header  is 
determined  from  the  following  expression: 


=  *^ej^VH  =  2.09  goDyuFyjjWtI  (24) 


where  the  factor  is  defined  as  the  vapor  header  occlusion  factor  for 

radiant  emission  to  space  and  given  a  value  of  0.85  for  this  analysis  (ref. 

13).  The  analysis  and  results  of  reference  13  (fig.  4)  although  primarily  for 
tubes  and  fins  are  assumed  valid  for  the  occlusion  of  the  header-panel  arrange¬ 
ment  shown  in  figure  1. 


Radiator  Tubes 

The  pressure  drop  in  the  radiator  tubes  where  flowing  vapor  is  condensing 
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TO-s  computed  from  a  combination  of  several  basic  flo-w  and  energy  equations  that 
are  given  in  detail  in  reference  6.  The  flow  model  (ref.  14)  used  assumed  that 
at  any  given  section  perpendicular  to  the  flow  direction^  the  temperature  and 
pressure  in  both  the  liquid  and  vapor  are  uniform  and  the  same  for  both  phases. 
This  flow  model  also  assumed  turbulent  flow  with  liquid  and  vapor  velocities 
uniform  in  each  phase  at  a  given  cross  section^  but  that  the  two  velocities 
were  not  necessarily  equal.  Pressure  drops  were  computed  for  a  series  of  in¬ 
cremental  tube  lengths^  and  the  pressure  drop  for  the  whole  tube  was  obtained 
by  summing  the  incremental  drops. 

The  total  change  in  pressure  for  an  entire  radiator  tube  is  comprised  of  a 
frictional  and  a  momentum  component.  The  friction  pressure  drop  is  described 
by  the  expression 


dPp  - 


-»|Re 


1.8 

g 


/o-092  tig 

\  S^fpg 


(25) 


where  Re  =  ^g  is  a  function  of  X.  The  differential  form  of 

the  change  in  pressure  due  to  a  change  in  momentum  is 


A  third  relation  is  required  that  relates  the  Increment  of  tube  length  and 
the  increment  of  condensate  formed: 


dW.= 


3600  ^  Jhp_  ^ 


7“^ 

h-r^  11.5  - 

Jg 


0.9 


dX 


+  _M  .  0-9^fe£^tot\ 

dWj. 


where  for  the  double  fin-tube  radiator  the  definition  of  K  is 

^b'HR 


K  =  4 


rtDi  eJ 


(27) 


(28) 


Equations  (25),  (26),  and  (27)  along  with  equation  (l)  are  solved  simultaneous- 
ly  for  dW^  ^^m^  total  change  in  pressure  for  each  increment  can 
then  be  foimd  from  the  relation 


dP  =  dP^  +  dPp 


(29) 


and  the  total  pressinre  change  for  the  entire  radiator  tube  can  be  found  by 
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summing  the  incremental  changes.  For  comparative  purposes^  the  total  tube 
pressure  drop  is  expressed  as  the  ratio  of  pressure  drop  to  the  tube  inlet 
pressure. 

The  pressiore  drop  associated  with  turning  the  vapor  from  the  vapor  header 
into  the  radiator  tubes  and  accelerating  the  flov  can  be  calculated  from  the 
expression 


entrance 


Ktt  h  p, 


where  the  tube  entrance  loss  factor 
tion  (l). 


is  given  the  same  value  used  in  eq^ua* 


Liq[uid  Header 

The  pressure  drop  in  the  liquid  header  is  obtained  by  applying  Fanning’s 
equation  with  a  friction  factor  for  turbulent  flow.  This  expression  is 

Q. 00102 


•where  Re  is  the  Reynolds  n-umher  corresponding  to  the  maximum  liquid  velocity 
Vrjj  that,  in  these  calculations,  -was  taken  as  4  feet  per  second.  The  liquid 
header  diameter  -which  is  assumed  constant,  is  determined  by  applying  the 

continuity  equation  at  the  header  exit  to  give  the  expression 


RADIATOR  WEIGHT  AND  GEOMETRY 

The  heat-rejection  analysis  of  the  fin  and  the  tube  'was  nondimensionalized 
so  that  the  results  could  be  used  for  general  design  purposes.  Using  the  pre¬ 
vious  results  for  analyzing  the  merits  of  the  variable  tube  -wall  double  fin- 
tube  configuration  in  radiator  designs  is  impractical  'without  consideration  of 
the  total  weight  of  the  vapor  and  liquid  headers  as  well  as  the  tubes  and  fins. 
It  is  necessary,  therefore,  to  consider  the  ratio  of  heat  rejection  per  unit 
weight  Ojpe-A  influencing  effects  of  radiator  geometry  limitations  on 

the  maximum  heat  rejection  per  unit  weight. 


Armor  Thickness 

In  order  to  determine  both  radiator  weight  and  geometry,  the  effects  of 
meteoroid  penetration  must  be  considered  on  both  the  tubes  and  headers.  This 
will  dictate  the  required  armor  protection  thickness  needed  for  the  radiator. 
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The  tuhe  armor  thickness  6^^  is  determined  using  the  meteoroid  protection  cri¬ 
teria  given  in  reference  9^  which  is  hased  on  a  comprehensive  appraisal  of  the 
available  data  and  theories  concerning  the  meteoroid  penetration  phenomenon. 
According  to  reference  9^  the  resultant  eq^uation  for  the  armor  thickness  & 
is  given  by  the  expression  ^ 


where  in  the  previous  equation 
a  =1.75 


P  =1.34 


pp  =  0.44  g/cu  cm 

Vp  =  98^400  ft/sec 

a  =  O.SSxlO"^^  g^/(sq  ft) (day) 

Insertion  of  these  constants  into  equation  (33)  along  with  utilizing  Young *s 
modulus  in  the  definition  of  sonic  velocity  of  the  material  yields  the  more 
compact . form 


6 


a 


1.48 


r  1 

-In  P(0) 


0.249 


(34) 


The  total  exposed  area  to  be  protected  by  direct  impacts  is  assmed 

to  be  the  outer  surface  of  the  vapor  header  and  the  projected  area  of  the  tube 
block.  The  liquid  header  contribution  is  assumed  to  be  negligible  since  its 
sinrface  area  is  small  compared  to  that  of  the  vapor  header.  Thus 

=  A(.  +  Ayjj  (35) 

The  radiator  tube  projected  area  is  given  by  the  expression 

\  -  (36) 

*  ^)’1S 

where  Z/R^  is  obtained  from  equation  (18)  and  Tjg  from  equation  (20).  The 

fraction  of  the  total  radiator  heat  that  is  rejected  by  the  tubes  and  fins 
is  defined  by  the  expression 

a  +  (X. 

^f  =  =  1  -  (37) 
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The  vulnerable  area  of  the  vapor  header  is  assumed  to  he  its  full  surface 
area  and  is  given  as 

a™  .  ^  (38) 

“'Vs 

•where  Fx^  is  the  occlusion  factor  for  the  effect  of  the  radiator  panel  on  the 
Vn 

heat  rejection  of  the  vapor  header  surface. 


Combining  equations  (35)  to  (38)  yields 


1  -  X. 


Weight  Ratio 

The  total  heat  rejection  per  unit  -weight  of  a  fin-tube  radiator  can  be  ex¬ 
pressed  as 

^re  j 

f40^ 

NZ 

where  for  radiation  from  both  sides  of  the  radiator  and  from  the  vapor  header 


2.09 


-  ae  4RQT^Tigll  +  r  )  + 


Insertion  of  eepiation  (23)  into  equation  (4l)  yields  the  following  form  of  the 
total  radiator  heat  rejection  per  xmit  length  of  tube: 


Qre.i  _  l\( .  *„4  ,  %H^VH^3 


where  the  outside  tube  diameter  (block  width)  can  be  expressed  in  terms  of 

the  liner  thickness  6^,,  the  inside  tube  diameter  and  the  armor  thickness 

8„  and  can  be  given  as 
a 

D  =  D.  +  28  +  25  (43 

o  1  c  a 


The  total  weight  of  the  radiator  is  comprised  of  the  individual  weights  of 
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the  vapor  header,  the  liquid  header,  and  the  tube-fin  panel.  The  weight  of  the 
liquid  inventory  in  the  subcooler  portion  of  the  radiator  tubes  is  neglected 
since  the  subcooler  length  is  small.  The  vapor  header  weight  is  given  by 
(ref.  6) 


W 


“  3  ^[(%H  +  ^VhJ^^Pc  +  (%H  +  25^^  + 


,)^aPaj 


and  that  of  the  liq^uid  header  and  condensate  hy 


"^LH  “ 


PjJ^LH 


Pc^H^(^H  +  ^LHp)  +  PaSa(% 


a\^E  +  ^a  +  26lh 


(44) 


(45) 


The  fin  and  tube  panel  weight  can  he  calculated  using  the  expression 


^tf  =  4ZtpjNZ  +  Nz[(Di  +  26^)2  -  pf]  +  NZp^[4R^R^  "  i  (^i  +  25^)^]  (46) 

When  t  =  al2T^/k]l^,,  Dq  from  equation  (43),  I  from  equation  (15),  and 

%  =  Rq  +  +_(Ss/Sa)^a  introduced  into  equation  (46),  the  equation  for 

the  tube  and  fin  panel  weight  becomes 


W. 


tf 


4p^NZaT^ 


+  6, 


+ 


8a)  2  +  PcJt6cWZ(Dq  +  6^) 


+ 

“D 

2(D.  +  26^  +  28^) 

%  -  (y 

rt(R. 


Sc) 


) 


(47) 


The  total  radiator  weight  is  then  obtained  by  summing  the  results  of  equa¬ 
tions  (44),  (45),  and  (47).  The  denominator  of  equation  (40)  can  be  found  by 
dividing  the  total  radiator  weight  by  the  total  tube  length  IJZ.  This  result 
along  with  the  results  of  equation  (42)  when  inserted  into  equation  (40)  yield 
the  radiator  heat  rejection  per  unit  weight.  The  peak  value  of  Q^e'/W  can 

then  be  obtained  by  plotting  the  resilLts  of  equation  (40)  as  a  function  of 
L/Ro  specific  values  of  conductance  parameter  N„  and  tube  side-wall  ra- 

tlo  S^/6^. 


Panel  Geometry 

In  addition  to  the  important  aspect  of  minimizing  weight  for  practical  ra¬ 
diator  designs,  it  is  also  of  Interest  in  most  cases  to  investigate  the  geom¬ 
etry  of  the  radiator  as  it  might  affect  the  integrating  of  the  vehicle  and  ra¬ 
diator.  Planform  area,  aspect  ratio,  and  fin  thickness  are  three  facets  of  the 
geometry  of  the  radiator  panel  that  must  be  determined  in  order  to  satisfy 
radiator -space  vehicle  Integration  and  structural  and  fabricational  require- 
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merits  of  the  fin-tuhe  configuration. 

Radiator  planform  area  Ap  is  obtained  from  the  equation 


A  =  2NZ(L  +  R^)  = 


(48) 


It  is  seen  from  this  equation  that  planform  area  will  vary  inversely  with  over¬ 
all  fin-tube  thermal  effectiveness  T]p  for  a  specific  choice  of  power  and  tem¬ 
perature  level.  The  planfoimi  area  will  generally  increase  with  increasing 
L/Rq  because  decreases  as  the  L/Rq  ratio  is  increased  (fig.  8). 

Another  interesting  factor  with  respect  to  the  geometry  of  the  radiator  is 
the  magnitude  of  the  fin  thickness.  Radiator  applications  might  require  that 
the  fin  have  structural  or  fabricational  qualities  that  could  result  in  non¬ 
optimum  weights  and  dimensions.  Fin  thickness  t  for  the  closed  sandwich  con¬ 
figuration  can  be  calculated  from  the  expression 


t 


oT^Z^ 

"ksy 


(49) 


The  panel  aspect  ratio ^  which  is  defined  as  the  ratio  of  panel  width  W 
to  tube  length  is  obtained  by  using  equation  (23)  for  W  and  the  results 
of  the  pressure  drop  calculations  for  Z. 


Method  of  Solution 

Simultaneous  solution  of  equations  (4)^  (20)^  (21)^  (29)^  (33)^  and  (39) 
requires  inputs  of  inside  tube  diameter  tube  and  vapor  header  pressure 

drop  ratios^  power  level,  temperature  level,  and  the  properties  of  the  mate¬ 
rials  and  cycle  fluid.  This  results  in  values  of  the  parameters  Tp, 

R,  and  uq  ^’or  a  selected  variation  of  L/Rq,  conductance  parameter  R^,  and 
the  tube  side-wall  ratio 

Important  results  required  for  weight  and  geometry  calculations  will  in¬ 
clude  the  number  of  tubes  R,  tube  length  Z,  tube  outside  radius  R^,  panel 
width  W,  and  the  inside  diameter  of  the  vapor  header  Dy^.  An  electronic 
digital  computer  was  used  to  obtain  the  desired  results  that  required  an  itera¬ 
tive  type  solution. 


RESULTS  ARD  LISCUSSIOR 
Calculation  Inputs 

Calculations  making  use  of  the  resultant  equations  developed  in  the  analy¬ 
sis  for  the  double  fin-tube  geometry  require  inputs  such  as  inside  tube  diam-  . 
eter,  radiator  vapor  inlet  temperature,  cycle  power  level  and  conditions,  mate- 
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rials  of  construction,  meteoroid  protection  criteria,  tube  block  side-wall 
thickness  ratio,  and  pressure  drop  in  the  tubes  and  header  in  order  to  com¬ 
pletely  specify  a  radiator  solution.  For  this  reason  two  specific  cases  are 
used  in  order  to  show  the  effects  of  some  of  the  previously  mentioned  variables 
on  heat  rejection  per  unit  weight  and  on  radiator  geometry. 

The  first  case  considered  is  a  1-megawatt  electrical  output  powerplant 
with  the  radiator  at  1700°  R.  The  tube  armor  and  the  fin  were  assumed  to  be 
made  of  beryllium,  and  the  tube  liner  was  assumed  to  be  a  columbium  alloy.  The 
second  case  is  a  500-kilowatt  system  with  the  radiator  at  1700°  R.  For  this 
case,  the  tube  armor,  tube  liner,  and  the  fin  were  all  taken  to  be  columbium  - 
1-percent  zirconium  alloy.  These  two  cases  both  used  a  peak  cycle  temperature 
of  2460°  R  and  potassium  as  the  cycle  fluid.  Inside  tube  diameters  of  3/8, 

1/2,  5/8,  3/4,  and  1  inch  were  chosen  with  tube  lengths  increased  to  allow 
100°  R  of  subcooling.  A  500-day  mission  time  and  a  probability  of  no  puncture 
P(0)  of  0.995  were  chosen  for  the  calculation  of  meteoroid  protection  thick¬ 
ness.  Pressure  drop  ratios  for  the  two  cases  were  set  at  AP/P  =  0.02  for  the 
vapor  header  and  AP/P  =0.05  for  the  radiator  tubes.  The  emittance  of  the 
surface  coating  on  the  fins,  tubes,  and  headers  was  taken  to  be  0.90.  Radiator 
material  properties  were  assumed  constant  with  temperature  and  evaluated  at 
1700°  R.  The  all  columbium  alloy  radiator  had  a  density  of  530  pounds  per  cu¬ 
bic  foot,  thermal  conductivity  of  34  Btu  per  hour  per  foot  per  °F,  and  a  modu¬ 
lus  of  elasticity  of  0. 202xl0l0  pounds  per  square  foot.  The  beryllium  radiator 
used  a  density  of  115  pounds  per  cubic  foot,  a  thermal  conductivity  of  51.5  Btu 
per  hour  per  foot  per  °F,  and  a  modulus  of  elasticity  of  0.397xl0l0  pounds  per 
square  foot.  The  beryllium  radiator  used  a  columbium  liner  with  the  material 
properties  previously  mentioned. 

Calculations  using  the  given  inputs  and  specifications  were  also  made  for 
the  central  fin-tube  geometry  that  used  the  method  described  in  reference  6 . 


Radiator  Weight 

Radiator  heat  rejection  per  unit  weight  C^g^/W  was  plotted  for  each 
value  of  inside  tube  diameter  Dp  and  tube  block'^  side -wall  ratio 
chosen  for  the  comparison  over  a  range  of  the  parameters  and  l/r  .  Re¬ 

sults  showing  the  variation  in  heat  rejection  rate  per  unit  weight  as  a  func¬ 
tion  of  l/Rq  for  several  values  of  conductance  parameter  are  shown 

plotted  for  two  sample  cases  in  figure  9.  The  figure  shows  results  for  the  tube 
block  side-wall  thickness  to  tube  armor  thickness  ratio  5s/&a  equal  to  0.5 
with  tube  inside  diameters  of  3/8  and  5/8  inch  chosen  for  the  500-kilowatt  and 
1-megawatt  cases,  respectively.  Each  constant  Ng  curve  is  seen  to  peak  at  a 

specific  value  of  L/Rq  with  the  value  of  the  L/Rq  at  peak  Qrej/^  increas¬ 
ing  as  Ng  is  increased. 

Plotting  the  maxima  for  each  Rg  curve  of  figure  9  and  the  additional  re¬ 
sults  for  the  Sg/Sg  =  0  and  1.0  cases  yield  a  performance  map  (fig.  lO)  that 
plots  peak  heat  rejection  per  unit  weight  against  L/r^.  It  is  seen  from  fig¬ 
ure  10  that  decreasing  the  armor  block  side-wall  thickness  resuAts  in  a  sub- 
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stantial  increase  in  the  value  of  Q^ej/W.  According  to  figure  10,  the  maxi- 
rniiTn  Qrej/^  occurs  at  an  Nq  approximately  equal  to  0.75  for  the  500-kilowatt 
system  at  a  tube  Inside  diameter  of  s/s  inch.  The  results  for  the  1-megawatt 
system  indicate  this  maximLim  also  occurs  near  an  Nq  =  0.  75  for  a  near  optimum 
tube  inside  diameter  of  5/8  inch. 

Plotting  the  maxima  of  figures  10(a)  and  (b)  on  figures  ll(a)  and  (b),  re¬ 
spectively,  and  the  results  for  additional  inside  tube  diameters  yield  the 
maximum  value  of  Qrej/^  corresponding  inside  tube  diameter  for  the 

three  values  of  6s/Sa  chosen  for  the  comparison.  The  curves  given  in  figure 
11(a)  for  the  500-kilowatt  case  indicate  that  maximum  QrejA  choice 

of  bg/^a  occurs  at  a  tube  inside  diameter  between  s/S  and  l/s  inch  with  the 
Qrej/w  curve  relatively  flat  in  this  region.  Maximum  Qrej/^  occ\irred  at  a 
diameter  between  l/2  and  5/8  inch  for  the  1-megawatt  system  as  shown  in  fig¬ 
ure  11(b). 

The  weight  results  of  the  double  fin-tube  configurations  are  compared  to 
the  central  fin-tube  results  in  figures  ll(a)  and  (b)  in  order  to  show  the  de¬ 
sirability  of  the  double  fin-tube  with  reduced  tube  side-wall  thickness.  It  is 
observed  from  figure  11  that  the  central  fin-tube  configuration  has  a  larger 
heat  rejection  per  unit  weight  than  the  double  fin-tube  with  bg/ba  =  1.0.  For 
the  500-kllowatt  case,  the  double  fin  tube  with  bg/Sa  =0.5  affords  an 
11-percent  weight  advantage  at  maximum  Qrej/^  over  the  central  fin  tube,  and 
a  maximum  upper  limit  of  39  percent  when  bg/Sa  ~  0*  radiator  weight  per 

kilowatt  of  electrical  power  for  the  500-kilowatt  case  using  a  columbium  radi¬ 
ator  is  10.2,  15.0,  and  19.0  pounds  per  kilowatt  for  the  bg/ba  =  0,  0.5,  and 
1.0  cases,  respectively. 

The  percentage  weight  savings  are  reduced  for  the  1-megawatt  power  level 
with  the  beryllium  radiator  and  are  8  percent  for  the  bg/ba  =0.5  case  and 
32  percent  for  the  bg/ba  =  case.  For  this  power  level,  the  specific  weight 
was  3.50,  4.68,  and  5.73  pounds  per  kilowatt  for  the  bg/ba  =  0,  0.5,  and  1.0 
cases,  respectively.  The  beryllium  radiator  has  a  smaller  percent  increase  in 
Qj.ej/W  than  the  columbium  radiator  because  there  is  less  relative  weight  in 
the  beryllium  armor. 

The  conclusions  reached  in  figure  19  of  reference  7,  which  compared  the 
double  fin-tube  and  central  fin-tube  configurations  without  taking  into  account 
the  effects  of  header  heat  rejection  and  weight  as  well  as  pressure  drops  in 
the  radiator  tubes,  showed  a  similar  percentage  increase  in  QrejA  reduc¬ 
ing  bg/bg^.  This  data,  which  was  for  a  1-megawatt,  1700°  R,  beryllium  radia¬ 
tor,  Indicated  a  6-percent  weight  advantage  for  the  double  fin  tube  at  a 
6s/ba  =0.5  and  a  28  percent  advantage  at  bg/bg  =  0  over  the  central  fin- 
tube  configuration. 

The  conductance  parameter  obtained  at  maximum  heat  rejection  per  unit 
weight  is  plotted  as  a  function  of  Inside  tube  diameter  for  a  bg/ba  =0.5  for 
both  geometries  in  figure  12.  Results  obtained  for  bg/bg  other  than  0.5 
showed  very  little  variation  in  the  optimum  value  of  Kc*  Both  the  500- 
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kilo-watt  and  1-mega-watt  examples  are  shown  with  the  double  fin-tuhe  results 
yielding  the  higher  value  of  for  a  specific  choice  of  inside  tube  diam¬ 
eter.  It  is  also  noted  that  the  values  of  obtained  for  the  500-kilowatt 

columbium  radiator  are  somewhat  larger  than  those  obtained  for  the  1-megawatt 
beryllium  radiator.  Reference  7  indicated  the  double  fin-tube  configuration 
reached  maximum  much  larger  values  of  than  those  obtained  in 

this  report. 


Radiator  Geometry 

Planform  area.  -  The  planform  area  of  the  double  fin-tube  geometry  is 
sho-wn  plotted  in  figure  13  for  the  two  power  levels  chosen  for  this  investiga¬ 
tion.  This  figure  illustrates  the  calculated  variations  of  planform  area  with 
L/Rq  ratio conductance  parameter  and  tube  side-wall  thickness  to  armor 

thickness  ratio  power  cycles  for  peak  Qrej/^  conditions. 

The  two  inside  tube  diameters  used  for  the  results  shown  in  figure  13  corre¬ 
spond  to  near  minimum  weight  conditions.  Reducing  the  tube  side-wall  thickness 
ratio  for  peak  L/Rq  ratios  results  in  only  a  small  variation  in  planform  area 
for  a  constant  value  of  but  results  in  a  sizable  decrease  in  L/Rq.  How¬ 

ever^  it  should  be  kept  in  mind  that  the  actual  fin  length  I  is  not  decreas¬ 
ing  that  rapidly  since  the  fin  length  is  being  increased  by  the  amount  6g^  -  bg 
as  6s/&a  approaches  zero.  The  decrease  in  L/Rq  (which  is  a  measure  of  the 
tube  center-to-center  distance)  as  bs/^a  decreases  causes  a  reduction  in  tube 
block  effectiveness  (due  to  reduced  oj^)  and  an  increase  in  fin  effectiveness 
(fig.  6)  that  result  in  only  a  small  variation  in  total  fin-tube  effectiveness 
for  a  specific  choice  of  This  accounts  for  the  negligible  variation  of 

planform  area  with  varying  bg/bg^.  Similar  results  are  obtained  for  the  other 
inside  tube  diameters  investigated. 

Comparison  of  the  planform  area  results  of  the  double  fin-tube  geometry 
with  those  of  the  central  fin-tube  geometry  indicates  good  agreement  as  shown 
by  figure  14.  The  planform  area  obtained  for  the  central  fin  tube  is  less  than 
that  of  the  double  fin  tube  for  internal  tube  diameters  greater  than  1/2  inch 
for  the  500-kilowatt  columbi-um  radiator.  At  the  3/8-inch  diameter^  which  cor¬ 
responds  to  near  minimum  weighty  the  two  fin-tube  geometries  agreed  to  within 
4  percent.  For  the  1-megawatt  beryllium  system,  the  double  fin  tube  has  the 
larger  planform  area  at  tube  diameters  greater  than  5/8  inch.  At  near  optimum 
weight  corresponding  to  a  5/8-inGh  diameter,  the  central  fin-tube  geometry 
offers  only  a  small  savings  in  planform  area.  The  planform  area  obtained  for 
the  double  fin-tube  geometry  in  reference  7  for  the  1-megawatt  radiator  is  much 
greater  than  that  obtained  from  the  comparison  calculations  of  this  report. 

This  is  brought  about  because  the  conductance  parameter  at  maximum  Qrej/^ 

the  case  of  reference  7  was  larger  than  the  values  obtained  in  this  report. 
Figure  13  verifies  that  increasing  conductance  parameter  Rq  results  in  in¬ 
creased  planform  area. 

Fin  thickness.  -  The  effect  of  reducing  the  tube  side -wall  ratio  bg/6g 
on  fin  thickness  is  shown  for  the  two  examples  in  figure  15  for  peak 
conditions  at  near  optimum  inside  tube  diameter.  In  both  cases,  reducing 
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6  /&  results  in  a  sizable  decrease  in  fin  thickness  at  constant  Accom¬ 

panying  this  decrease  in  fin  thickness  is  a  reduction  in  L/Rq  and  thus  tuhe 
spacing.  The  magnitude  of  the  fin  thickness  obtained  for  the  two  sample  cases 
at  optimum  values  of  \  (fin  thickness  t  is  greater  than  0.020  in.  for  the 
columbium  radiator  and  greater  than  0.035  in.  for  the  beryllium  radiator)  re¬ 
sulted  in  values  of  fin  thickness  that  are  reasonable  and  should  satisfy  struc¬ 
tural  and  fabricational  requirements. 

The  comparison  of  the  fin  thickness  obtained  for  the  double  fin-tube  con¬ 
figuration  indicates  an  increasing  total  fin  thickness  with  increasing ^values 
of  tube  side  wall  to  armor  ratio  Bg/Sa  throughout  the  range  of  investi¬ 

gated  for  both  power  levels  (fig.  16).  Curves  are  given  for  fin  thickness  at 
maximum  -/^  both  geometries  and  both  system  power  levels  with  central 

and  double  fin-tube  Bg/Bg^  =  0.5  cases  agreeing  closely  for  the  1  megawatt 
radiator . 

Panel  aspect  ratio.  -  Radiator  panel  aspect  ratio,  which  is  defined  as  the 

ratio  of  panel  width  W  to  the  tube  length  Z  for  the  four  panel  radiator  of 
figure  1,  is  shown  plotted  in  figure  17  for  the  two  double  fin-tube  example 
cases.  The  aspect  ratio  of  a  panel  is  seen  to  decrease  as  the  tube  block  side- 
wall  ratio  Bg/Sg  decreases  at  a  constant  Ng.  The  decrease  in  aspect  ratio 
with  reduced  Sg/6a  primarily  a  result  of  the  reduced  L/Rq  since  reduc¬ 
tions  in  Bg/Bg  had  little  effect  on  planform  area  as  shown  by  figxare  13.  ^ 
These  results  can  be  indicated  by  inspection  of  equation  (48),  which  describes 
planform  area  and  length,  and  equation  (23),  which  describes  the  total  panel 
width  W.  The  reduction  in  aspect  ratio  at  constant  planform  area  is  explained 
by  the  calculations  which  showed  that  as  Bg/Sg,  gets  smaller  the  number  of 
tubes  and  the  tube  length  increase  as  the  L/Rq  ratio  decreases.  This  was 
accompanied  by  a  large  reduction  in  the  working  fluid  inlet  tube  velocity. 

The  panel  aspect  ratio  showed  practically  no  difference  between  the  cen¬ 
tral  fin-tube  and  double  fin-tube  configurations  as  indicated  in  figures  18(a) 
and  (b)  for  the  500-kllowatt  and  1-megawatt  systems,  respectively. 

Number  of  tubes.  -  The  number  of  radiator  tubes  required  was  found  to  in¬ 
crease  substantially  as  the  tube  block  side-wall  ratio  decreased.  The  500- 
kilowatt  columbium  radiator  with  a  tube  inside  diameter  of  3/8  inch  had  298, 
337,  and  430  tubes  for  Sg/Bg  ratios  of  1.0,  0.5,  and  0,  respectively,  at 
maximum  heat  rejection  per  unit  weight.  The  1-megawatt  beryllium  radiator 
using  a  tube  inside  diameter  of  5/8  inch  had  223,  250,  and  302  tubes  for  ^a 

ratios  of  1.0,  0.5,  and  0,  respectively,  at  maximum  Q^ej/’^- 

For  simplicity  and  reliability  of  fabrication  it  is  desirable  to  reduce 
the  number  of  tubes.  A  sizable  reduction  in  the  number  of  tubes  without  much 
reduction  in  '^gj/w  can  be  had  by  just  increasing  the  size  of  the  tube  inside 
diameter.  For  the  500-kilowatt  columbium  radiator  with  tube  block  side-wall 
ratio  set  at  0.5,  a  3/8-inch  inside  tube  diameter  required  337  tubes  whereas  a 
l/2-lnch  inside  diameter  results  in  206  tubes.  This  reduction  in  the  number  of 
tubes  results  in  only  a  2-percent  decrease  in  1-megawatt  beryl- 
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lium  radiator  with  5g/&^  set  at  0.5  required  250  tubes  for  a  5/8-inch  inside 
tube  diameter  whereas  a  3/4-inch  inside  diameter  required  just  182  tubes.  Ihis 
reduction  resulted  in  only  a  3-percent  decrease  in  ^g./W.  Similar  reductions 

®  nmber  of  tubes  can  be  obtained  at  values  of  6  /&  equal  to  0  and  1.0 
vith  small  reductions  in  Q  ./W.  ^  ^ 

_  Comparison  of  the  number  of  tubes  for  the  central  and  double  fin-tube 
radiators  at  maximum  i^gj/w  for  the  500-kilowatt  case  indicated  the  central 
fin-tube  geometry  required  308  tubes  compared  to  298  tubes  for  the  double  fin- 
ube  Og/Og^  -  1.0  radiator  at  an  inside  tube  diameter  of  3/8  inch.  This  trend 
also  held  for  the  1-megawatt  case  for  which  the  central  fin-tube  required  226 
tubes  compared  to  223  for  the  double  fin-tube  8g/8g  =1.0  radiator. 

of  results.  -  If  the  values  of  the  ratio  of  5  /&  other 
than  the  three  values  chosen  for  this  investigation  (O,  0.5,  and  1)  txe  desig¬ 
nated  by  forthcoming  meteoroid  bumper  results,  additional  points  may  be  ob- 
ained  from  curves  faired  through  the  calculated  points  to  obtain  intermediate 
es  of  panel  aspect  ratio  w/z,  fin  thickness  t,  radiator  planforra  area 
Ap,  and  radiator  heat  rejection  per  unit  weight  Qrej/W- 


SUMMAEY  OF  RESULTS 

An  analysis  of  the  double  fin-tube  configuration  with  variable  tube  block 
SI  e-wall  thickness  and  a  comparison  with  a  comparable  central  fin-tube  radia¬ 
tor  for  the  sample  cases  of  a  1-megawatt  beryllium  radiator  and  a  500-kilowatt 
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A  substantial  increase  in  radiator  heat  rejection  per  unit  weight  is 
possible  for  the  double  fin-tube  geometry  if  the  tube  block  side-wall  thickness 
can  be  reduced  to  less  than  75  percent  of  that  required  for  normal  armor  thick- 
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j-v  1^"  'double  fin-tube  configuration  with  ratio  of  tube  block  side-wall 
thickness  to  normal  armor  thickness  of  0.5  offers  an  11  and  8  percent  weight 
savings  over  the  central  fin-tube  geometry  for  the  500-kilowatt  and  1-megawatt 
cases,  respectively.  This  corresponds  to  15.0  pounds  per  kilowatt  for  the  500. 
kilowatt  case  and  4.68  pounds  per  kilowatt  for  the  1-megawatt  case  considered 

pounds  per  kilowatt  for  the  central  fin-tube 
ouu-kilowatt  and  1-megawatt  cases,  respectively. 


The  percent  weight  savings  is  Increased  to  39  and  32  percent  for  the 
500-kilowatt_and  1-megawatt  cases,  respectively,  if  the  ratio  of  tube  block 
side-wall  thickness  to  normal  armor  thickness  is  set  at  0.  In  this  case  the 
corresponding  radiator  specific  weights  would  be  10.2  and  3.50  pounds  per  kilo¬ 
watt  for  the  double  fin-tube  geometry  and  16.7  and  4.83  pounds  per  kilowatt  for 
the  central  fin-tuhe  geometry. 


4.  The  double  fin-tube  configuration  reached  maximum  heat  rejection  per 
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unit  weight  for  the  1-megawatt  system  at  a  conductance  parameter  of  0.75  where¬ 
as  maximum  heat  rejection  per  unit  weight  for  the  central  fin-tuhe  geometry 
occurred  at  0.63.  These  values  for  the  500-kilowatt  case  were  0.63  and  0.75, 
respectively,  for  the  double  and  central  fin-tuhe.  However,  the  conductance 
parameter  associated  with  maximum  heat  rejected  per  unit  weight  decreases  with 
decreasing  inside  tube  diameter  for  both  the  double  and  central  fin-tube  radi¬ 
ators. 

5.  The  double  fin-tube  radiator  reaches  maximum  heat  rejected  per  unit 
weight  for  a  tube  Inside  diameter  in  the  range  from  3/8  to  l/2  inch  for  the 
500-kllowatt  case,  and  from  l/2  to  5/8  inch  for  the  1-megawatt  case. 

6.  The  physical  dimensions  of  the  double  fin-tube  radiator  (planform  area, 
panel  aspect  ratio,  number  of  tubes,  inside  tube  diameter,  and  fin  thickness) 
can  be  varied  over  a  fairly  wide  range  without  seriously  decreasing  the  radia¬ 
tor  heat  rejection  per  unit  weight. 

7.  The  fin  thickness  obtained  for  maximum  heat  rejection  per  unit  weight 
for  the  double  fin-tube  geometry  is  of  reasonable  fabricational  and  structural 
magnitude  (greater  than  0.02  in.  for  the  500  kw  columbium  radiator  and  greater 
than  0. 035  in.  for  the  1  Mw  beryllium  radiator) . 

8.  The  values  of  planform  area  and  panel  aspect  ratio  obtained  for  the 
double  and  central  fin-tube  geometries  were  in  close  agreement  regardless  of 
the  choice  of  tube  block  side-wall  ratio  near  the  minimum  weight  condition. 
Comparison  of  the  total  fin  thickness  for  the  two  geometries  agreed  well  at  a 
tube  block  side-wall  thickness  ratio  of  0.5,  but  showed  considerable  varia¬ 
tions  at  side-wall  ratios  of  1.0,  and  0. 

9.  Reduction  of  the  tube  block  side-wall  thickness  ratio  substantially 
increases  the  number  of  radiator  tubes  at  maximum  heat  rejection  per  unit 
weight.  However,  the  nimiber  of  tubes  can  be  reduced  to  the  values  correspond¬ 
ing  to  the  central  fin-tube  geometry  by  increasing  the  value  of  above 

that  for  maximum  heat  rejection  per  unit  weight.  The  weight  penalty  involved 
is  only  several  percent . 


Lewis  Research  Center 

National  Aeronautics  and  Space  Administration 
Cleveland,  Ohio,  August  12,  1964 
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Figure  4.  -  Assumed  radiator  heat- transfer  paths. 
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Figure  5.  -  Fin  temperature  profile;  closed  sandwich  fin  tube 


Ratio  of  fin  half-length  to  tube  outer  radius,  L/R^ 

(b)  Conductance  parameter,  1.0  and  10. 

Figure  6.  -  Pin  effectiveness;  double  fin-tube  geometry.  Powerplant 
output,  500  kilowatts;  temperature  at  radiator  inlet,  1700°  R- 
Inside  tube  diameter,  3/4  inch;  columblum  radiator.  ' 
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Figure  8.  -  Overall  fin-tube  effectiveness;  double  fin- tube 

Figure  7.  -  Tube  effectiveness;  double  fin- tube  geometry.  Powerplant  output,  500  kilowatts;  temperature 

geometry.  Powerplant  output,  500  kilowatts;  at  radiator  inlet,  1700®  R;  inside  tube  diameter,  3/4 

temperature  at  radiator  inlet,  1700°  R;  in-  inch;-  columbium  radiator, 

side  diameter,  3/4  inch;  columbium  radiator. 


Maximum  heat  rejection  per  unit  weight 
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(a)  Powerplant  output,  500 
kilowatts;  columblum 
radiator. 


(b)  Powerplant  output,  1  megawatt;  beryl 
lliom  radiator. 


Figure  11.  -  Maximum  heat  rejection  per  unit  weight.  Temperature  at  radiator 
inlet,  1700  R 
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Radiator  fin  thickness  at  peak  heat  rejection  per  unit  weight;  double  fin-tube  geometry.  Temperature  at  radiator  inlet 

1700^  R. 
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(a)  Powerplant  output,  500  kilowatts;  Inside 
columblum  radiator. 
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Figure  17.  -  Radiator 
unit  weight;  double 
at  radiator  Inlet, 


panel  aspect  ratio  at  peak  heat 
fin-tube  geometry;  four  panels. 
1700°  R. 
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**T!he  aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knoivU 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof.” 

— ^National  Aeronautics  and  Space  Act  of  1958 
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